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Abstract 

This is a brief survey of laws of large numbers, fluctuation results and 
large deviation principles for asymmetric interacting particle systems that 
represent moving interfaces on the plane. We discuss the exclusion pro- 
cess, the Hammersley process and the related last-passage growth models. 



1 Introduction 

This article is a brief overview of recent results for a class of stochastic processes 
that represent growth or motion of an interface in two-dimensional Euclidean 
space. The models discussed have in a sense rather orderly evolutions, and the 
word "directed" is included in the title to evoke this feature. 

Let us begin with generalities about these stochastic processes. The state 
at time t G [0, oo) is of the form h{t) = {hi{t) : i e Z) with the interpretation 
that the integer- or real- valued random variable hi{t) represents the height of 
the interface over site i of the substrate Z. We call the state h = {hi) a height 
function on Z. The interface on the plane is then represented by the graph 
{(«, hi) : i G Z}. Each particular process has a state space that defines the set of 
admissible height functions. The state space will be defined by putting restric- 
tions on the increments (discrete derivatives) hi — of the height functions. 

The random dynamics of the state are specified by the jump rates of the 
individual height variables hi. The rates are functions of the current state h. 
For the sake of illustration, suppose that only jumps ±1 are permitted for each 
variable hi. Then the model is defined by giving two functions p{h) and q{h). 
If the current state is h, then the height value ho at the origin jumps down to a 
new value h'o = ho — 1 with rate p{h), and jumps up to a new value h^ = hg + 1 
with rate q{h). In the spatially homogeneous case the rates for hi are p{6ih) 
and q{9ih) where the spatial translations are defined by {Oih)j = hi^j. 

The quantity p{h) is the rate for the jump ho r\ ho — 1 in an instantaneous 
sense: in the current state h, the probability that this jump happens in the 
next infinitesimal time interval (0, dt) is p{h)dt A- o{dt). Rigorous constructions 
of the processes utilize Poisson processes or "Poisson clocks." A rate A Poisson 
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process N{t) is a simple continuous time Markov chain: it starts at A''(0) = 0, 
runs through the integers 0, 1,2,3,... in increasing order, and waits for a rate A 
exponential random time between jumps. A rate A exponential random time is 
defined by its density ip{t) = Xe~^* on R+. The number of jumps N{s+t) — N{s) 
in time interval (s, s + has the mean Xt Poisson distribution 



P{N{s + t)- N{s) =k} = 



k\ 



{k > 0). 



If the overall rates are bounded, say by p{h) < A, then a Poisson clock with a 
time-varying rate p(h(t)) can be obtained from N{t) by randomly accepting a 
jump at time t with probability p{h(t))/X. 

Later we mention in passing rigorous constructions of some processes. In 
each case the outcome of the construction is that all the random variables 
{hi{t) : i € Z, i > 0} are defined as measurable functions on an underlying 
probability space P). Since these processes evolve through jumps, the 

appropriate path regularity is that with probability 1, the path t i— > h{t) is 
right-continuous with left limits (cadlag for short). The use of Poisson clocks 
makes the stochastic process h{t) a Markov process. This means that if the 
present state h{t) is known, the future evolution (/i(s) : s > t) is statistically 
independent of the past {h{s) : < s < t). This is a consequence of the "for- 
getfulness property" of the exponential distribution. For a complete discussion 
of these foundational matters we must refer to textbooks on probability theory 
and stochastic processes. 

This article covers only asymmetric systems. Asymmetry in this context 
means that the height variables hi{t) on average tend to move more in one 
direction than the other. For definiteness, we define the models so that the 
downward direction is the preferred one. In fact, the great majority of the 
paper is concerned with totally asymmetric systems for which q{h) = 0, in 
other words only downward jumps are permitted. Symmetric systems behave 
quite differently from asymmetric systems, hence restricting treatment to one 
or the other is natural. 

Stochastic processes with a large number of interacting components such 
as the height process h{t) = {hi{t) : i € Z ) belong in an area of probability 
theory called interacting particle systems. (jSpitzeii . Il970l) is one of the semi- 
nal p apers of this subject. Here is a selection of books and lecture notes on th e 
topic: 



De Masi and Presuttil (ll99lh.lDurretlJ Jlissh.lKipnis and LandimI d 19991) . 

■Liggett! (ll985h.lLiggettl (li999h . lLiggettl (|2004l ). ]Varadhanl (|200(i r" Krug and SpohnI 
(ll992^ and'Spohnl (jl99ll ) are sources that combine mathematics and the theo- 



retical physics side. 

Our treatment is organized around three basic questions posed about stochas- 
tic models: (i) laws of large numbers, (ii) fluctuations, and (iii) large deviations. 

(i) Laws of large numbers give deterministic limit shapes and evolutions 
under appropriate space and time scaling. A parameter n oo gives the ratio 
of macroscopic and microscopic scales. A sequence of processes h"(t) indexed 
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by n is considered. Under appropriate hypotheses the height process satisfies 
this type of result: for {t,x) € M+ x M 

n~^h"^^-^{nt) u{x,t) as n ^ oo, (1) 

and the limit function u satisfies a Hamilton- Jacobi equation ut + ,f{ux) = 0. 

(n) Fluctuations. After a law of large numbers the next question concerns 
the random fluctuations around the large scale behavior. One seeks an exponent 
a that describes the magnitude of these fluctuations, and a precise description 
of them in the limit. A typical statement would be 

h? , (nt) — nulx, t) 
'""'^ ' -— — . Zit, x) (2) 

where Z{t, x) is a random variable whose distribution would be described as 
part of the result. The convergence is of a weak type, where it is the probability 
distribution of the random variable on the left that converges. 

(hi) Large deviations. The vanishing probabilities of atypical behavior fall 
under this rubric. Often these probabilities decay as e~^"'^ to leading order, 
with another exponent p > 0. The precise constant C G (0, oo) is also of interest 
and comes in the form of a rate function. When all the ingredients are in place 
the result is called a large deviation principle (LDP). An LDP from the law of 
large numbers ([T]) with rate function / : R ^ [0, oo] could take this form: 

lim \im n-HogP{h?,{nt)e{z-e,z + e)} = -I{z) (3) 

valid for points z in some range. Positive values I{z) > represent atypical 
behavior, while limit ^ would force I{u{x,t)) — 0. 

Example. For classical examples of these statements let us consider one- 
dimensional nearest-neighbor random walk. Fix a parameter < p < 1. Let 
{Xfe} be independent, identically distributed (IID) ±l-valued random variables 
with common distribution P{Xk = 1} = p = 1 — P{Xk = —1}. Define the 
random walk by Sq = 0, Sn = Sn-i + Xn for n > I. Then the strong law of 
large numbers gives the long term velocity: 

lim n^^Sn — V where v — E(Xi) = 2p — 1. 

n—*oc 

The convergence in the limit above is almost sure (a.s.), that is, almost every- 
where (a.e.) on the underlying probability space of the variables {^fc}. 

The order of nontrivial fluctuations around the limit is n^^^ ("diffusive") and 
in the limit these fluctuations are Gaussian. That is the content of the central 
limit theorem: 



.^f Sn — nv 1 , , , 1 f 



e 

oo 



-^'/'dz. 
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The parameter = E[(Xi — u)^] is the variance. 
Random walk satisfies this LDP: 

Um hm - logPIIS*,! - nx\ < ne} = -I(x) (4) 

where the rate function /: M ^ [0,oo] is identicaUy oo outside [—1, 1] and 

^ 1-2;, 1 — x 1 + a; , 1 + x ^ r -, n 

/(x) = ^log^^ + ^log^ for. e [-1,1]. (5) 

I{x) can be interpreted as an entropy. Convex analysis plays a major role 
in large deviation theory. Part of the general theory behind this simple case 
is that / is the convex dual of the logarithmic moment generating function 
A(6') = logE(e^^i). 

Results for random wa lk are cover ed in graduate probability texts such as 
( Durretdliool and (K allenbergL l2002h . 



At the outset we delineated the class of models discussed. Important models 
left out include diffusion limited aggregation (DLA) and first-passage percola- 
tion. Their interfaces are considerably more complicated than interfaces de- 
scribed by height functions. But even for the models discussed our treatment is 
not a complete representation of the mathematical progress of the past decade. 
In particular, this article does not delve into the recent work on Tracy- Widom 
fluctuations, Airy processes and determinantal point processes. These topics are 
covered by many authors, and we give a number of references to the literature 
in Sections 13.11 and 13.21 Overall, the best hope for this article is that it might 
inspire the reader to look further into the references. 

Recurrent notation. The set of nonnegative integers is Z+ = {0, 1,2,...}, 
while N = {1, 2, 3, . . . }. The integer part of a real x is [x] = maxjn E Z : n < x}. 
a\/ b — max{a, b} and a Ah = min{a, b}. 

Acknowledgement. The author thanks M. Balazs for valuable consultation 
during the preparation of this article. 



2 Limit shape and evolution 

We begin with the much studied corner growtli model and a description that 
is not directly in terms of height variables. Attach nonnegative weights {^ijj} 
to the points of the positive quadrant N'^ of 1?, as in Figure [TJ Yij 

represents the time it takes to occupy point (?, j) after the points to its left 
and below have been occupied. Assume that everything outside the positive 
quadrant is occupied at the outset so the process can start. Once occupied, a 
point remains occupied. Thus this is a totally asymmetric growth model, for 
the growing cluster never loses points, only adds them. 

Let G{k,£) denote the time when point (k,£) becomes occupied. The above 
explanation is summarized by these rules: G{k,i) = for {k,£) ^ N^, and 

G{kJ)^G{k-l,e)VG{k,£-l)+Yk^i ioi{k,£)eN^. (6) 
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Figure 1: Each point G N has a weight Yij attached to it. 



The last equality can be iterated until the corner (1, 1) is reached, resulting in 
this last-passage formula for G: 



G{k,l) 



max 



(7) 



(ij")67r 



where Tik,i is the collection of nearest-neighbor up-right paths tt from (1, 1) to 
(fc,£). Figure [21 represents one such path for (fc,£) = (5,4). 

This model and others of its kind are called directed last-passage percolation 
models. "Directed" refers to the restrictions on admissible paths, and "last- 
passage" to the feature that the occupation time G{k, £) is determined by the 
slowest path to {k,£). (By contrast, in first-passage percolation occupation 
times are determined by quickest paths.) 



1 2 3 4 5 
Figure 2: An admissible path from (1, 1) to (5,4) 
Our first goal is to argue the existence of a limit for n^^G{[nx], [ny]) as 
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n — > oo. Assume now that the weights {^ij} are IID nonnegative random 
variables. 

The idea is to exploit sub (super) additivity. Generalize the definition of 
G{kJ) to 

G{{k,£),{m,n)) = max K^j, 

(ij)e7r 

where Tl(^k,e).{m,n) is the collection of nearest-neighbor up-right paths tt from 
(k + 1,£ + 1) to (m, n). The definitions lead to the superadditivity 

G{k, e) + G{{k, £), (m, n)) < G{m, n). (8) 



Kingman's subadditive ergodic theorem (jPurrettl . |2004 Chapter 6) and some 



estimation implies the existence of a deterministic limit 

'y{x,y) = lim n^^G{[nx], [ny]) for {x,y) G M^. (9) 



Moment assumptions under whi ch the limit fun ction 7 is continuous up to 
the boundary were investigated bv iMartin ( 2004 ). We turn to the problem 



of computing 7 explicitly, and for this we need very specialized assumptions. 
Essentially only one distribution can be currently handled: the exponential, 
and its discrete counterpart, the geometric. Take the {i^ij} to be IID rate 1 
exponential random variables. In other words their common density is e~^. 

The difficulty with finding the explicit limit has to do with the superaddi- 
tivity. The limit in Birkhoff's ergodic theorem is simpl y the expectati on of the 
function averaged over shifts: X^Li f ° ^ E/ dPurrettl . 12004 Chapter 
6). But the subadditive ergodic theorem gives only an asymptotic expression 
for the limit. We need a new ingredient. We shall embed the last-passage model 
into the totally asymmetric simple exclusion process (TASEP). This has explic- 
itly identifiable invariant distributions ("steady states") with which we can do 
explicit calculations. 

Originally TASEP was introduced as a particle model. We wish to link 
TASEP with the last-passage model in a way that preserves the original formu- 
lation of TASEP, while mapping particle occupation variables into height incre- 
ments and particle current into column growth. To achieve this we transform 
the coordinates of Figure[T]via the bijection (i,j) 1-^ (« — j, — j)- The result 
is the last-passage model of Figure [31 Weights are relabeled as Xij = Yi-j-j. 
The transformation of admissible paths is illustrated by Figure [H Let the new 
last-passage times be denoted by H{k,£). For £ < A fc the maximizing-path 
formulation uses now paths of the kind represented in Figure |31 

H{k,e) = max V X,j (10) 

where E/c^^ is the collection of paths cr from (0, —1) to (fc,£) that take steps of 
two types: (1,0) and (—1,-1). The connection with the previous last-passage 
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Figure 3: Last-passage model for TASEP. The horizontal i-axis and the vertical 
j-axis are labeled, and points on the i-axis from —3 to 3 are marked. Weights 
Xij are attached to points («, j) such that i E Z, j E — N and j < A i. 

process is H{k,£) = G{k — £, —£). The process {H{k,£)} is also defined by the 
recursion 

H{kj.) = H{k^i,e)y H{k + i,e + i) + Xk^i {£<OAk) (ii) 

together with the boundary values H{k,£) = for ^ > A fc. It satisfies the 
limit 



X{x,y) = limn ^H{[nx],[ny]) ~ j{x — y,—y) for ?/ < A a;. (12) 



Figure 4: The image of the path in Figure[21 Now it goes from (0, —1) to (1, —4). 



To establish the TASEP connection we first define a height process w{t) = 
{wi{t) : i E Z) that will turn out to be an alternative description of the last- 
passage process {H{i,j)}. Initially at time t ~ the height is given by 




(13) 
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This is the boundary of the region {j < A i} filled with Xij^s in Fig. [31 This 
initial shape is a wedge, hence the symbol w. 

Give each column i an independent rate 1 Poisson clock Ni. Variable Wi 
jumps downward according to this rule: if i is a jump time for Poisson process 
Ni, then 

/ N , ^ . , , I Wi-i{t—) = Wi(t—) — 1 and ,^ 

w,{t) = w,{t-) - 1 pYOvided { ' (14) 

Equivalently, each column variable Wi jumps down independently at rate 1, as 
long as the state w{t) remains in the state space 

Xi^ {hell- -.hi - e {0, 1} for all i £ Z}. (15) 

(See Fig. [5] for an example.) The interaction between the variables is encoded 
in rule p^ . It forces the timc-cvolution of each variable Wi to depend on the 
evolution of its neighbors. By contrast, without the interaction the variables Wi 
would simply march along as Poisson processes independently of each other. 



j 
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Figure 5: A possible height function w{t) (thickset graph) with column values 
W-i{t) = —3, wo{t) = —2, wi(t) = —2, etc. 8 jumps have taken place during 
time (0,t]. The columns grow downward, x's mark the allowable jumps from 
this state. 

This construction defines the height process w{t) = {'Wi{t) : i e Z) for all 
times t £ [0, cx)) in terms of the family of Poisson clocks {Ni}. Given this process 
w(t) define the stopping times 

T{i,j) = M{t>0:w,it) <j} (16) 

that mark the time when column i first reaches level j. Stopping time is a 
technical term for a random time whose arrival can be verified without looking 
into the future. 

Initial condition (flSl) imphes T{i, j) — for j > i AO. Rule (fT4)) tells us that 
T(i — 1, j) V T{i + l,j + 1) is the stopping time at which the system is ready 
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for Wi to jump from level j + 1 to j. (Note that Wi must have reached level 
j + 1 already earher because if the rules are followed, T{i + 1) comes after 
T{i,j + 1).) By the forgetfulness property of the exponential distribution, after 
the stopping time T{i —l,j)\/T{i + l,j+l) it takes another independent rate 1 
exponential time Xij until Wi jumps from j + 1 to j. Consequently the process 
{T{i,j)} satisfies the recursion 

T{t,j) = T{i - V T{t + + X,,,. (17) 

This is of the same form as the recursion ([TT]) satisfied by {H{i,j)}. From 
this one can prove that indeed the processes {T(i, j)} and {H{i,j)} are equal 
in distribution. Therefore gives also n~^T{[nx],[ny]) X{x,y). This is 
the precise meaning of the earlier claim that the height process w{t) gives an 
alternative description of the last-passage process {H{i,j)}. 

Subadditivity and some estimation justifies the existence of a concave func- 
tion g on R such that 

lim t~^wixt]{t) = g{x) a.s. for x G M. (18) 

t — >oo 

Since rates are 1, g records only the initial height outside the interval [—1,1], 
and so 

g[x) = A X for \x\ > 1. 

Since the interface is a level curve of passage times, X(x, g{x)) = 1 for —1 < x < 
1. The last-passage limits are homogeneous in the sense that X{cx, cy) = cX{x, y) 
for c > 0. Consequently A and then 7 can be obtained from g. 

To summarize, thus far we have converted the original task of computing 
j{x,y) of ^ to finding the function g of (HH]) on the interval [—1,1]. Now 
consider the general height process h{t) with state space Xi from pS]) and 
dynamics defined as for w{t) above: height variable hi jumps one step down at 
every jump epoch of the Poisson clock Ni, provided this jump does not take the 
height function out of ^"1. A jump attempt that would violate the state space 
restriction is simply ignored. 

Certain technical issues may trouble the reader. An infinite family of Poisson 
clocks has infinitely many jumps in any nonempty time interval (0, e). So there 
is no first jump attempt in the system and it is not obvious that the local 
rule leads to a well-defined global evolution: to determine the evolution of hi 
on [0,i] we need to look at the evolution of its neighbors /ii±i, and then their 
neighbors, ad infinitum. However, given any T < 00, almost surely there are 
indices \ —00 and ij, y 00 such that iV^^ and TVjj^ have no jumps during 
(0,T]. Consequently the system decomposes into (random) finite pieces that do 
not communicate before time T. The evolution can be determined separately in 
each finite segment which do experience only finitely many jumps up to time T 
(again almost surely). Another technical point is that the clocks {Ni} have no 
simultaneous jumps (almost surely) so one never needs to consider more than 
one jump at a time. 
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Given that the height process h{t) has been constructed, next define the 
increment process r](t) = {rji{t) : i € Z) by 

r/,(i) = h,{t) - h,.i{t). (19) 

Process 'q{t) has compact state space {0, 1}^ and its dynamics inherited from 
h{t) can be succinctly stated as follows: each 10 pair becomes a 01 pair at 
rate 1, independently of the rest of the system. To sec this connection, observe 
that if Poisson clock Ni jumps at time t, the height process undergoes the 
transformation hi{t) — hi{t—) — 1 only if {r]i{t—) = l,77i+i(i— ) = 0), and 
then after the jump the situation is {rii{t) — 0,7/i+i(t) — 1). This is a direct 
translation of the condition that jumps are executed only if the state h remains 
in the state space Xi. 

It is natural to interpret the I's as particles and the O's as holes, or va- 
cant sites. The process ri{t) is the totally asymmetric simple exclusion process 
(TASEP). In this model the only interaction between the particles is the exclu- 
sion rule that stipulates that particles are not allowed to jump onto occupied 
sites. This property is enforced by the evolution because the definitions made 
above ensure that a jump in Poisson clock Ni sends a particle from site i to 
site i + I only if site i + 1 is vacant. Total asymmetry refers to the property 
that particles jump only to the right, never left. The definitions also entail this 
connection between the heights and the particles: 

hi{0) — hi{t) ~ cumulative particle current across the edge (i, i + 1). (20) 

We need to discuss two more properties of these processes, (i) stationary 
behavior and (ii) the envelope property. Then we are ready to compute the 
function g of (fT5|) . 

(i) Stationary behavior. For p G [0, 1], the Bernoulli probability measure Vp 
on {0, 1}^ is defined by the requirement that 

Upitj ■.n, = l for iel,rij = {) for j G J} = p'^' (1 - p)''^' (21) 

for any disjoint I, J with cardinalities |/| and | J|. Measure corresponds 
to putting a particle at each site independently with probability p. 

It is known that the measures {i^p}pe[o,i] are invariant for the process r]{t), 
and in fact they are the extremal members of the compact, convex set of invari- 
ant probability measures that are also invariant under spatial shifts. Invariance 
means that if the process rjit) is started with a random i/p-distributed initial 
state 77(0), then at each time t >0 the state r]{t) is t-p-distributed, and further- 
more, the probability distribution of the entire process rj{-) ^ {ri{t) : t > 0) is 
invariant under time shifts. 

If wc know the current state h{t), then the probability that jumps down in 
a short time interval {t,t + e) is er]o{t){l — r]i{t)) + 0{e^). This follows because a 
jump can happen only when a 10 pair is present, and from properties of Poisson 
processes. Estimation of this kind proves that 

ho{t) - hoiO) = - f ryo(t)(l - mit))ds + M{t) (22) 
"'0 
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where M{t) is a mean-zero martingale. This identity is a stochastic "fundamen- 
tal theorem of calculus" of sorts. Since things are random the difference between 
ho{t) — ho{0) and the integral of the infinitesimal rate cannot be identically zero. 
Instead it is a martingale. This is a process whose increments have mean zero 
in a very strong sense, namely even when conditioned on the entire past. 

Let us average over ((^ in the stationary situation. Let Ej^^ denote expecta- 
tion of functions of the stationary process r]{-) whose state ri{t) is i/p-distributed 
at each time t. Normalize the height process h{-) at time zero so that /io(0) = 0. 
Then /i(-) is entirely determined by r]{-). Since 77o(s) and 771 (s) are independent 
at any fixed time s, we get 

E,^[ho{t)]^^tf{p) (23) 

where the particle Rux is defined by 

/(p) = p(l-p). (24) 

(a) Envelope property. Even though the flux / is nonlinear and therefore, as 
we see later, TASEP is governed by a nonlinear PDE, the height process has a 
valuable additivity property. Suppose a given initial height function h{0) € Xi 
is the envelope of a countable collection {z''-^\G)} keK of height functions in the 
sense that 

/i,(0) = sup zf ^ (0) for each site i e Z. (25) 

Take a single collection {Ni\ of Poisson clocks, and let all processes h{t)^ z'^^\t) 
evolve from their initial height functions by following the same clocks {Ni\. 
This kind of simultaneous construction of many random objects for the purpose 
of comparison is called a coupling. By induction on jumps one can prove that 
this coupling preserves the envelope property for all time: 

Lemma 2.1 hi{t) — sup zl''\t) for all i eT, and all t>Q, almost surely. 

We take as auxiliary processes z'^^^t) suitable translations of the basic wedge 
process defined in p^ - p^ . For fc £ Z set 

^f'(0) = |'~^' and z«(t)^/.,(0)-K-f (t). (26) 

10, i> k 

The apex of the wedge z'^^\Q) is at the point (fc, hk{0)), and then the definition 
of the wedge ensures that /i(0) > z('^^(0). Hypothesis (|^ holds and Lemma [TT] 
gives this variational equality: 

h,{t)^snY>{hu{Q)+wf\t)}. (27) 

fcez 

Now we extract two results from the assembled ingredients: the function 
g of (fT8|) . and a general "hydrodynamic limit" that describes the large scale 
evolution of the process. 
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First specialize (P7|) to the stationary situation where ho{0) = and the 
increments are i/p-distributed, and write (|27|) in the form 

t-'ho{t) = sup{i-i/i[,,](0) + t-ii«([*^l)(i)}. (28) 

yeR 

Let t — > oo. Inside the braces on the right (0) —> py a.s. by the law 

of large numbers. t~^w'^^^'^\t) 9{^y) by a translation of the limit (fT8|) . 
With some work take the limit outside the supremum. Then we know t~^ho{t) 
converges. By supplying some moment bounds we can take expectations over 
the limits, and with (|23p arrive at 

-/(p) =sup{py + g(-y)}. (29) 

yew 

This is a convex duality (jRockafellaii . Il970l) . From the explicit invariant distri- 
butions (|2ip we obtained / in and then we can solve (P^)) for g. (Without 
the invariant distributions we can carry out part of this reasoning but we cannot 
find / and g explicitly.) Let us record the results. 

Theorem 2.2 For the limit ^ g{x) = -^(1 - x)^ for -1 < a; < 1. For the 
limit © ^{x,y) = [./x + for x,y >Q. 

We turn to the hydrodynamic limit. Assume given a function ito on M and 
a sequence of random initial height functions /i"(0) ^ Xi {n E N) such that 

n~"^/i"j^](0) uo{x) a.s. as n ^ oo for each a; g R. (30) 

For this to be possible uq has to be Lipschitz with < UQ(a;) < 1 Lebesgue-a.e. 

Theorem 2.3 For a; S R and t > we have the limit 

n~^h"j^^-^{nt) u{t,x) a.s. as n —s- oo (31) 

where 

u{t,x) = sup{uo{y) + tg(^—^)}. (32) 

yeR ^ \ t / J 



Equation ([5^ is a Hopf-Lax formula ( Evand . 1998 ) and it says that 
entropy solution of the Hamilton-Jacobi equation 



ut + /(ux)=0, u\t= 



Uq. 



is the 



(33) 



In other words this equation governs the macroscopic evolution of the height 
process. Theorem 12. 31 is proved by showing that, as n — + oo, variat ional formula 
P7|) for n~^h? Ant) turns into (15^ . Details can be found in ( SeppalainenI , 
1993). 
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F urthe r remarks. The function g in Theorem 12.21 was first calculated by 



Rostj (jl98lh in one of the seminal papers of hydrodynamic limits, but without the 
last-passage representation and with a different approach than the one presented 
here. 

Let us discuss various avenues of generalization. We encounter immediately 
difficult open problems. 

(i) Generalizations that retain the envelope property. The argument sketched 
above that combines the envelope property with the duality of th e flux and the 
wedge s hape to derive hydrodynamic limits was introduced in (|SepDalainenl . 



1998allc . 1999). An earlier instance of the variational connection appeared in 



Aldous and Diaconid ( 19951 ) for Hammer s ley's process. This work itself was 
based on the classic paper (jHammerslevl . 119721) : see Sect ion 13.21 below. Also 
in qu eueing literature similar variational expressions arise (jSzczotka and Kellv , 
19901) . 

To define the K-exclusion process we replace the state space Xi of ([15]) with 



Xk ^{heZ^ ■.0<h,- h,^i < K for all i e Z} 



(34) 



for some 2 < ii' < oo. Otherwise keep the model the same: rate 1 Poisson clocks 
{Ni\ govern the jumps of height variables /i^, and jumps that take the state h 
outside the space are prohibited. The increment process is now called to- 
tally asymmetric iiT-exclusion (some authors use "generalized exclusion"). The 
variational coupling (Lemma 12. ip works as before. But invariant distributions 
are unknown, and there is even no proof of existence of an extremal invariant 
distribution for each density value p G [0, if]. No alternative way to compute 
/ and g has been found. Theorem 12.31 is valid, but the most that can be said 
about / and g is that they exist as concave functions. 

Interestingly, the situation becomes again explicitly analyzable for K = oo 
where the only constraint on h is hi < /i^+i. The increment process is a spe- 
cial case of a zero range pr ocess. I t s stat e space i s (Z^ )^ and IID geometric 



distributions are invariant ( Liggettl . 119731 : lAndiell . Il982h . As a final step of 
generaliza tion, away from monotone height functions, let us mention bricklayer 



processes ([Balazsl l2003t iBalazs et al.l 120071 ) whose increments rji ^ hi — hi-i 



can be positive or negative. 

The variational coupling of Lemma 12.11 works equally well for certain mul- 
tidimensional height proce sses h(t) = (hi(t) : i g Z'^) of the type disc ussed 
here. Examples appear in (|Rezakhanloul . l2002bl : ISeppalainenL lioool l2007t ). No 
explicit invariant distributions are known for multidimensional height models. 
The variational scheme proves that scaled height processes converge to solutions 
of Hamilton- Jacobi equations as in Thcorem l2.3l But again one can only assert 
the existence of / and g instead of giving them explicitly. 

Another direction of generalization is to let the weights {i^i.j} have distri- 
butions other than exponential or geometric. The height process h{t) ceases 
to be Markovian but the last-passage model of Figs. [T] and [2] makes sense. As 
mentioned, the limit ^{x, y) in ([9]) is explicitly known only for the exponential 
and geometric cases. A distribution as simple as Bernoulli {Yi.j takes only two 
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values) cannot be handled. However, if the paths are altered to require that one 
or both coordinates increase strictly, th en the variat i onal approac h does find 
the exphcit shape for the Bernoulli case (jSeppalainenl Il997l Il998bl ). 

Thus the present situation is that an explicit limit shape can be found only 
for some fortuitous combinations of path geometries and weight distributions. 

(a) Partially asymmetric models. Let us next address the case where the 
column variables hi are allowed to jump both up and down. Fix two parameters 
< q < p such that p + q = 1 (convenient normalization). Give each column 



i two independent Poisson clocks, N, 



(-) 



with rate p and A^, 



(+) 



with rate q. At 

jump times of N^^' variable hi attempts to jump to hi ± 1, and as before, a 
jump is completed if its execution does not take the state out of the state space 
Xi. For the increment process this means that a 10 pair becomes a 01 pair at 
rate p, and the opposite move happens at rate q. Bernoulli distributions (j2ip 
are still invariant. This increment process is the asymmetric simple exclusion 
process (ASEP). In the same vein one can allow K particles per site and talk 
about asymmetric i^T-exclusion. 

T he envelope propert y of Lemma |2. II is now lost. An alternative approach 
from (jRezakhanloul . 120011 ) utilizes compactness of the random semigroups of the 
height process. Limit points are characterized as H amilton-Jacobi semigroups 
via the Lions-Nisio theorem ( Lions and Nisiol . 1982 ). Thereby Theorem 12.31 is 
derived for one-dimensional asymmetric if -exclusion. For K — 1 the flux / in 
must be replaced by f{p) = {p — q)p{^ — p), while for 2 < i^T < oo the 
flux is unknown. In the multidimensional case it is not known if the resulting 
equation itself is random or not. 



3 Fluctuations 

A simple way to create initial height functions /i"(0) that satisfy assumption 
(|30p is to take independent increments with distributions 

PK(0) = l]=«^(^)-wo(^)), 

and at the origin assign the deterministic value /ip (0) ~ ['^'"o(O)]. The stationary 
situation is of this type with uq{x) — px. Then initial fluctuations 

n-'^Hhl„,^iO) - nuoix)} 

are Gaussian in the limit n oo. This makes it natural to look for a distribu- 
tional limit at later times t > on the central limit scale n^^^: 

n^^/'^{h''^^^^{nt) - nu{x,t)} — > Cit,x) as n ^ cx), (35) 

for some limit process ({t,x). Such limits can be proved but process {C{t,x)} 
turns out to be a deterministic function of the initial fluctuations {C(0, x)}. Con- 
sequently limit p5|) docs not record any fluctuations created by the dynamics. 
Theorem 13.31 below gives a precise statement of this type. 
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In asymmetric systems the fluctuations created by the dynamics occur on 
a scale smaller than n^/^. Two types of such phenomena have been found. 
Processes related to the last-passage model and exclusion process discussed 
in Section [5] have order n^/^ fluctuations whose limits are distributions from 
random matrix theory. A class of linear processes has order n^/* fluctuations 
and Gaussian limits related to fractional Brownian motion with Hurst param- 
eter H = 1/4:. To see these lower order fluctuations one can start the system 
with a deterministic initial state, or one can start the system in the stationary 
distribution or some other random state but then follow the evolution along 
characteristic curves of the macroscopic PDE. The fluctua tion situation is very 
differ ent for symmetric systems; the reader can consult (jKipnis and Landiml . 
19991 Chapter 11). 



3.1 Exclusion process 

Probability distributions from random matrix theory were discovered as limit 
laws for last-passage growth models almost a decade ago. 

Theorem 3.1 ()johanssonl bOOOl ). For the corner growth model 

G{[xn], [ny]) - n-f{x,y) 



< s 



where F is the Tracy-Widom GUE distribution. 



Fis) 



oo, 



(36) 



The distribution F first appeared as the limit distributi on of the scaled 
large st eigenvalue of a random Hcrmitian matrix from the GUE ( Tracv and WidomI . 
19941 ) . GUE is short for Gaussian Unitary Ensemble. This means that a random 
Hermitian matrix is constructed by putting IID complex-valued Gaussian ran- 
dom variables above the diagonal, IID real-valued Gaussian random variables 
on the diagonal, and letting the Hermitian property determine the entries below 
the diagonal. Then as the matrix grows in size, the variances of the ent r ies ar e 
scaled appropriately to obtain limits. The standard reference is (|Mehtal . [200l . 

Theorem 13.11 and related results initially arose entirely outside probability 
theory (except for the statements themselves), involving the RSK correspon- 
dence and Gessel's identity from combinatorics and techniques from integrable 
systems to analyze the asymptotics of the resulting determinants. The RSK 
correspondence, named after Robinson, Schensted and Knuth, is a bijective 
mapping between certain arrays of integers or integer matrices (in this case the 
matrix in Figure [T] if the Yij^s are integers) and pairs of Young tableau x. These 
latter objects are u biquitous in combinatorics. Standard references are (jFulton . 
1997 ; Saganl . 200 ih . More recently determinantal point processes have appeared 
as the link between the growth processes and random matrix theory. We shall 
not pursue these to p ics further for many exce l lent r evie ws are ayailabl e: iBaik 
( 2005^ ■ iDeiftI (|2000^ ■ Ijohanssoni (|2002h . iKonid (|2005^ and lSpohnI (l2006h . 

Precise limits such as (|36l) have so far been restricted to totally asymmetric 
systems. Next we discuss ideas that fall short of exact limits but do give the 
correct order of the variance of the height for partially asymmetric systems. 
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Consider the height process h{t) whose increraents r]i{t) — hiit) ~ hi^i{t) 
form the asymmetric simple exclusion process (ASEP). This process was intro- 
duced in the remarks at the end of Section O Each height variable hi attempts 
downward jumps with rate p and upward jumps with rate g, and p > q. A 
jump is suppressed if it would lead to a violation of the restrictions ft-i-i < hi 
or hi > hi^i — 1 encoded in the state space Xi of (fT5|) . In the increment process 
each 10 pair becomes a 01 pair at rate p and each 01 pair becomes 10 pair at 
rate q. 

On large space and time scales the height process obeys the Hamilton-Jacobi 
equation 

Ut + f{ux) = with f{p) = (p - q)p{l - p). 

This PDE carries information along the curves i — f'{ux{t,x)), in the sense 
that the slope Ux is constant along these curves as long as it is continuous. At 
constant slope — p the characteristic speed is V — f'{p) = {p— — 2p). 

Consider the stationary process: < p < 1 is fixed, and at each time t > Q 
the increments {?7i(i)}igz have Bernoulli i^p-distribution from (|2ip . Normalize 
the heights by setting initially /io(0) — 0. We determine the order of magnitude 
of the variance of the height as seen by an observer traveling at speed V. 



Theorem 3.2 ( Balazs and Seppalainenl . l2007bl ) Height fluctuations along th 



e 



characteristic satisfy 

< liminf t~2/^Var{/iryp4i(t)} < limsupt"^/^Var{/i[vpfi(t)} < oo. 

t^oo 



If the observer choses any other speed v ^ V^, only a translation of initial 
Gaussian fluctuations would be observed. Take v > to be specific. Due to 
the normalization ho(0) = we can write 

On the right the first expression in parentheses is a height increment along a 
characteristic and so by the theorem has fluctuations of order t-^^^. The last sum 
of initial increments has Gaussian fluctuations of order t^^^ and consequently 
drowns out the first term. 

The proof of Theorem 13.21 is entirely different from the proofs of Theorem 
13.11 As it involves an important probabilistic idea let us discuss it briefly. 

Couplings enable us to study the evolution of discrepancies between pro- 
cesses. In exclusion processes these discrepancies are called second class parti- 
cles. Consider two initial ASEP configurations 77(0), C(0) £ {0, 1}^. The config- 
urations differ at the origin: ((0) has a particle at the origin (Co(0) = 1) but 
77(0) docs not (7?o(0) = 0). At all other sites i ^ we give the configurations 
a common but random value 77^(0) = Ci(0) according to the mean p Bernoulli 
distribution. Let the joint process {ri{t),({t) ■ t > 0) evolve together governed 
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by the same Poisson clocks. The effect of this coupling is that there is always 
exactly one site Q{t) such that CQ(t)(^) — 1; ''lQ(t){t) = 0, and r]i{t) — (i{t) for 
all i ^ Q{t). 

Q{t) is the location of a second class particle relative to the process ri{t). 
(Relative to ({t) one should say "second class antiparticle.") In addition to 
ordinary exclusion jumps, Q yields to ry-particles: if an yy-particle at Q + 1 
jumps left (rate q) then Q jumps right to switch places with the ?7-particle. 
Similarly an ry-particle at Q — 1 switches places with Q at rate p. These special 
jumps follow from considering the effects of clocks Nq2i on the discrepancy 
between 77 and 

Proof of Theorem 13.21 utilizes couplings of several processes with different 
initial conditions. Evolution of second class particles is directly related to dif- 
ferences in particle current (height) between processes. On the other hand Q 
and the height variance are related through this identity: 



Var{/i[,4](i)} - p{l - p)E{ \Q{t) - [vt] \ } for any 



(38) 



The right-hand side can be expected to have order smaller than t precisely when 
V — VP on account of this second identity: 



Eg(t) = tyP. (39) 

From these ingredients the bounds in Theorem 13.21 arise. 

Further remarks. As already suggested at the end of Section [21 a major 
problem for growth models is to find robust techniques that are not dependent 
on particular choices of probability distributions or path geometries. Progress on 
fluctuations of the corner growth model beyond the exponential case has come 
in situations that are in som e sense extreme: for distributions with heavy tails 
Hamblv and Martin . 20071) or for points close to the b oundary of the q uadrant 



Baik and Suidanl . 120051: iBodineau and Martini . l2005f ). See review bv iMartinI 



20061) . 



The second class particle appears in many places in interacting particle sys- 
tems. In the hydrodynamic limit picture second class particles converge to 
characteristics and shocks of the macr oscopic PDE ( iFerrari and Font"ei . 1994bl : 
Rezakhanloul Il995t ISeppalainenl . I2OOI ). Versi ons of identities (1551) and (IMl) ar e 
valid for zero range and bricklayer processes (jBalazs and Seppalainenl . l2007al ). 
Equation (|39l) is surprising because the process as seen by the second class 
particle is not stationary. 

In general the view of the process from the second class particle is compli- 
cated. Studies of invariant distributions seen by second class p articles appear 



( Derrida et all . Il993l: iFerrari et all . 11994 iFerrari and Martini . [ioOHl . There 



are special cases of parameter values for certain processes where unexpected 
simplification takes place and the process seen from the second class parti cle 
has a product-form invariant distribution ([Derrida et al.l 1X997 : lBalazsll200l . 
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3.2 Hammersley process 



We began this paper with the exclusion process because this process is by far the 
most studied among its kind. It behooves us to introduce also Hammersley's 
process for which several important results were proved first. It has an ele- 
gant graphical construction that is related to a classical combinatorial question, 
namely the maximal lengt h of an increasing s ubsequenc e of a random permu- 
tation. Th is goes back to ( Hammersley . 1972 ): see also ( Aldous and Diaconisl 
Il995lll999l) . 

We begin with the growth model. Put a homogeneous rate 1 Poisson point 
process on the plane. This is a random discrete subset of the plane characterized 
by the following property: the number of points in a Borel set B is Poisson 
distributed with mean given by the area of B and independent of the points 
outside B. Call a sequence {xi,ti), {x2,t2), ■ ■ ■ , {xk,tk) of these Poisson points 
increasing ii xi < X2 < ■ ■ ■ < Xk and ti < t2 < ■ ■ ■ < t^. Let L{{a,s), (&, t)) 
be the maximal number of points on an increasing sequence in the rectangle 
(a, 6] X (s, t] (Fig.H]). The random permutation comes from mapping the ordered 
x-coordinates to ordered ^-coordinates in the rectangle, and L((a, s), (&, t)) is 
precisely the maximal length of an increasing subsequence of this permutation. 
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Figure 6: Increasing sequences among planar Poisson points marked by x's. 
-L((a, s), (b, t)) = 4 as shown by the circled Poisson points that form an increasing 
sequence. 



The limit 



lim n^'^L{{0,0),{nx,nt)) = 2Vxt 



(40) 



holds for x,t > 0. The limit exists by superadditivity exactly as for ([9]). 
The functional form c^/xi follows from scalin g properties of the Poiss on pro- 
cess. The value c = 2 was first derived by IVersik and Kerov ( 19771 ) , while 
Logan and SheppI (1977) independently proved c > 2. The fluctuation result 
for L is ana l ogous to p6p . with normalization n^^^ and the Tracy- Widom limit 
(|Baik et al.Ul999f) . 
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We embed the increasing sequences in the graphical construction of the Ham- 
mersley process. This process consists of point particles that move on M by 
jumping. Put a rate 1 Poisson point process on the space-time plane and place 
the particles initially on the real axis. Move the real axis up at constant speed 
1. Each Poisson point {x,t) instantaneously pulls to x the next particle to the 
right of X. We label the particles from left to right: Zi{t) € K is the position 
of particle i at time t. We could regard the variables Zi as heights again, but 
the particle picture seems more compelling. This construction is illustrated by 
Fig. [T] In terms of infinitesimal rates, the construction realizes this rule: inde- 
pendently of other particles, at rate Zi — Zi-i variable Zi jumps to a uniformly 
chosen location in the interval {zi^i, Zi). 



time t 



Zi-l{t) Zi{t) 



Zi+l{t) Zi+2{t) 



time 



l(0) 



^^(0) 



Figure 7: Portion of the graphical construction of Hammersley's process, x 's 
mark space-time Poisson points, "'s mark particle locations at time and at a 
later time t > 0. Space-time trajectories of particles are shown. The horizontal 
segments are traversed instantaneously and the vertical segments at constant 
speed 1. 

As in Section [21 there is a variational characterization for this construction. 
Define an inverse for the maximal path variable L{{a, s), {b,t)) by 

r((a, s),t, w) = inf{/i > : L((a, s), (a -|- h, t)) > w}. 

Take an initial particle configuration {2^(0)} G that satisfies Zi-i{0) < Zi{0) 
and i~'^Zi{0) ^ as i — > — oo. Then the graphical construction leads to a 
well-defined evolution {zi{t)} that satisfies 

z,{t)^ inf {zfe(0)+r((zfc(0),0),t,z-fc)}. (41) 

k:k<i 

Here is the hydrodynamic limit. Consider a sequence of processes z"{t) that 
satisfies ri"-'^z^^^j(0) uo{y) for each y e M, say in probability. The initial 
function uq is nondecreasing, locally Lipschitz and satisfies 

y~^uo{y) -> as 2/ ^ -oo. (42) 
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Define ^ 

u{t,x) ^ inf {uo{y)+ |, (i, x) e (0, cx^) x R. (43) 

y:y<x L 4t J 

Since rates are unbounded, we need to assume a left tail bound to prevent the 
particles from disappearing to — oo: given e > there exist < g, no < oo such 
that 

P{z"(0) < —ei'^/n for some i < —qn} < e for n > uq. 
Under these assumptions 

n~^zjjj^] (ni) ^ u(i, x) in probability (44) 

( Seppalainen . 19961 ). The function defined by the Hopf-Lax formula solves 



the Hamilton- Jacobi equation ut + [ux)^ — 0. 

Let us state a precise result about the central limit scale fluctuations ([35]) 
that covers also shocks. For {t, x) S (0, oo) x M let 

I{t, x) = {ve (-00, x] : uit, x) = uoiy) + } (45) 



be the set of minimizers in (j43[) , guaranteed nonempty and compact by hypoth- 
esis (UH). Then {t,x) is a shock if I{t,x) is not a singleton. This is equivalent 
to the nonexistence of the x-derivative Ux{t,x). 

Fluctuations on the scale n^/^ from the limit ([33]) are described by the pro- 
cess 

Ut,x) = n-'/^{z;^,^^{nt) - nu{x,t)}. 

Assume the existence of a continuous random function Co on R such that the 
convergence in distribution Cn(0, • ) ^ Co holds in the topology of uniform con- 
vergence on compact sets. Deflne the process C by 

({t,x) = inf Co{y) 

y£l(t,x) 



where I{t,x) is the (deterministic) set defined in ([IS]) . 

Theorem 3.3 For each (t,x), (^nit,x) Citi^) distribution. 

As stated in the beginning of Section [S] this distributional limit reflects no 
contribution from dynamical fluctuations as the process C is a deterministic 
transformation of (o- The underlying reason is that the dynamical fluctuations 
of order n^/^ are not visible on the n^/^ scale. The dynamical fluctuations are 
the universal ones described by the Tracy- Widom laws. See again the discussion 
and references that follow Theorem 13. II 

Further remarks. The polynuclear growth model (PNG) is another re- 
lated (l-l-l)-dimensional growth model used by several authors for studies of 
Tracy- Widom fluctuations and the A i ry process i n the KPZ scahng picture 
(iBaik and Raind . I2OO0I: iFerraril 120041 ; jjohanssonl I2OO3 I; IPrahofer and Spohnl . 



2002ll2004^ ." Like the Hammersley process, the graphical construction of the 
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PNG utilizes a planar Poisson process, and in fact the same underlying last- 
passage model of increasing paths. This time the Poisson points mark space- 
time nucleation events from which new layers grow laterally at a fixed speed. 
Roughly speaking, this corresponds to putting the time axis at a 45 degree angle 
in Figs. [S] and [71 

More about the phenomena related to Theorem 13.31 can be f ound in article 
(ISeppalainen . 2002 ). A similar theorem for TASEP appears in ( Rezakhanloul. 



2002a ) . Earlier work on the diffusive fluctuations of ASEP was done by Ferrari and Fontei 
l|l994allbh . 



3.3 Linear models 

We turn to systems macroscopically governed by linear first order equations 
Ut + bux = 0. Fluctuations across the characteristic occur now on the scale n^^^ 
and converge to a Gaussian process related to fractional Brownian motion. 

Th e random average process (RAP) was first studied bv lFerrari and Fonted 
( 1998h . The state of the process is a height function cr : Z ^ R with (Ti G K. 
denoting the height over site i. (More generally the domain can be Z'^.) The 
basic step of the evolution is that a value ct^ is replaced by a weighted average of 
values in a neighborhood, and the randomness comes in the weights. This time 
we consider a discrete time process. The basic step is carried out simultaneously 
at all sites i. 

Now precise formulations. Let {it(fc,r) : fc S Z, t G N} be an IID collection 
of random probability vectors indexed by space-time Z x N. In terms of coor- 
dinates u{k,T) — {uj{k,T) : ~M < j < M). We assume the system has finite 
range defined by the fixed parameter M. We impose a minimal assumption that 
guarantees that the weight vectors are not entirely degenerate: 

P{maxMj(0,0) < 1} > 0. (46) 

3 

For technical convenience we also assume that the process is "on the correct 
lattice": there does not exist an integer h > 2 such that for some b € Z the 
mean weights p{j) = Eiij(0,0) satisfy J2j<Eb+hi.P(j) = 1- 

To start the dynamics let cr(0) be a given random or deterministic initial 
height function. The process <t(t), r = 0, 1, 2, . . . , is defined iteratively by 

<^iiT) =^Uj{i,T)a^+j{T^l), T>l,ieZ. (47) 

3 

As before, we can define the process of increments ?7i(T) — ai{T) — ai-i{T). 
The increments also evolve via random linear mappings and are conserved like 
particles in exclusion processes. 

As in Section[2]we create suitable initial conditions for a hydrodynamic limit. 
Consider a sequence of processes tr" (r) indexed by n e N, initially normalized 
by ctq (0) = 0. For each n assume independent initial increments {f^f (0) : i G Z} 
with 

E[r/f(0)] =p(i/n) and Yar[7]'^ (0)] ^ v{i/n) (48) 
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for given Holder 1/2 + e functions p and v. Assume a uniform moment bound: 
sup„., E[ |?7r(0)|2+'5 ] < (X) for some 5>Q. 

The hydrodynamic limit is rather trivial for it consists only of translation. 
Define a function u on M by u(0) = and u'{x) = p{x). The characteristic 
speed is 



Then for each [t, x) £ R+ x 



n 



"'^(Tr^^i ([nt]) — > u{x ~ ht) as n ^ c», in probability. 



In other words, the height obeys the linear PDE ut + hux = 0. 

On the central limit scale n^/^ one would also see only translation of initial 
fluctuations. To see something nontrivial we look at fluctuations around a char- 
acteristic line. Fix a point y G K. and consider the characteristic line t ^-^ y -'rtb 
emanating from (y,0). Define space-time process 

Zn{t.r) = ^[:,,]+[.^]+[„,,](H) - '^[:,,]+[.^](0) , {t,r) e M+ X M. 

The spatial variable r describes fluctuations around the characteristic on the 
spatial scale n^/^. For the increment process 0) represents net current 

from right to left across the characteristic. 

Theorem 3.4 ( Balazs et al. . 20061) The finite- dimensional distributions of the 
process n^^^^Zn converge to those of the Gaussian process {z{t,r) : t > 0, r € 
R} described below. 

The statement means that for any finite collection of space-time points 
(ti,ri), . . . , {tk,rk), the M'^-valued random vector n~i/**(Z„(ti, ri), . . . , Zn{tk,rk)) 
converges in distribution to the vector {z{ti,ri), . . . , z{tk,rk))- The limiting pro- 
cess z has the following representation in terms of stochastic integrals: 



z{t,r) = p{y)<JaVK (p^2(t„5)(r-z)dW^(s,z) 

+ V'viv) / sign(x-r)$^2t(-|a;-r|)dB(x). 



Above is a 2-parameter Brownian motion on M+ x M and S is a 1-parameter 
Brownian motion on M independent of W. The first integral represents dynam- 
ical noise generated by the random weights, and the second the initial noise 
propagated by the evolution. The functions in the integrals are Gaussian den- 
sities and distribution functions: 

1 ( x"^ } /■^ 

V<T2(a;) = ^==cxp<^ -— ^ ^ and $^2(0;)=/ ip„2{y)dy. 

The only effects from the initial height are the mean p{y) and variance v(jj) of the 
increments around the point ny. The parameter ct^ is the variance of the prob- 
abilities p{j), and K another parameter determined by the distribution of the 

weights. Process z has a self-similarity property: {z(ai, a^/^r)} = {a^/'^z{t,r)}. 
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In the special case where v{y) = K.p{yY the temporal process {z(i,r) : t G 
K+l (for any fixed r) has covariance 



Ez(s, r)z{t, 



{^s + ^^t- ^J\t- s\ 



This identifies z{- ,r) as fractional Brownian motion with Hurst parameter 
H — 1/4. In particular, this limit arises in a stationary case where the av- 



eraging involves two points and the weight is beta distributed (jBalazs et al 
liooa Example 2.1). 



The proof of Theorem 13.41 utilizes a special case of another stochastic model 
of great contemporary interest, namely random walk in random environment 
(RWRE). Here is how the RWRE arises. An environment uj = {u(fc,T)} is 
determined by the weight vectors. Given ui, define a "backward" walk {Xl''^ : 
s £ Z+} on Z with initial position Xq'^ = i and transition probability 

P'^iXl^;, =y\Xl^^^x)^ Uy^.,ix, r - s), s = 0, 1, 2, . . . 

The superscript oj on P'^ indicates that it is the quenched path measure of 
{Xl'"^ : s € Z_|_} under a fixed oj. The basic step (|47|) of RAP evolution can be 
rewritten so that (Ti^r) equals the average value of the previous height function 
(7(r — 1) seen by a walk started at i after one step: 



This can be iterated all the way down to the initial height function: 

a,(T) = i?-[a^.,.(0)]. 

Note that the expectation E'^ over paths of the walk XI' under fixed weights 
to sees the initial height function {(7^(0) }igz as a constant. 

We have here a special type of RWRE called "space-time." Another term 
used is "dynamical environment" because after each step the walk sees a new 
sample of its environment. Proof of Theorem 13.41 requires limits for t he walk it- 
self and its quenched mean process E'^ iXt.-' '^ ). These results appear in (jBalazs et al 
20061: iRassoul-Agha and Seppalainenll2005h . 



Independent walks on Z display the same behavior as RAP. Let the process 
Zn{t,r) be the net particle current across the characteristic t i— > [ny] + [r^/n] 
+ [ntb] where b is the common average speed of the particles. Then under suit- 
able assumptions on the initial particle arrangements and their jump kernel, Zn 
satisfies a stronger form of Theorem 13.41 that also contains process-level conver- 
gence. One adjustment is necessary: the constants in front of the stochastic 
integrals in are different for the random walk case. The stationary system 
sees again fractional Brownian motion as the fimit of the current . Details for 
the random walk case appear in ( Seppalainen . 20051 : Kumar . 2007 ). Earlier re- 
lat ed results for a P oisson system of independent Brownian motions appeared 
in (|Diirr et all . Il985f) . 
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4 Large deviations 



We present the large deviation picture for the Hammersley process, so we 
continue in the setting of Section 13.21 Recall the definition of the longest 
path model among planar Poisson points illustrated in Fig. [B) Abbreviate 
Ln = i((0, 0), (n, n)). Then the limit (gO]) is n^^Ln 2. Here is the large 
deviat ion theorem for L n- It was completed shortly before the fluctuation re- 



sult of iBaik et al 



Logan and Shepn 
(|1999^ . 



1998 ). through a c ombination of severa l ind ependent papers: 



1977 ). Kim ( 1996f ). Seppalainen ( 1998d) and Deuschel and Zeitouni 



Theorem 4.1 We have the following upper and lower tail large deviation bounds. 

lim n^^ logP{L„ > nx} = -I{x) for x > 2 (50) 

n — >oo 

with rate function I{x) — 2a::cosh~"'^(a;/2) — 2\/ x^ — 4 . 

lim n^"^ logP{L„ < nx\ = -U{x) for0<x<2 (51) 

n — *C30 

with rate function U{x) — R2{s) ds where i?2(s) — slog(s/2) — s + 2 is the 
rate function for IID mean 2 Poisson random variables. 

To develop this theme further we state a lower tail LDP for the tagged par- 
ticle in Hammersley's process. An interesting feature is that the large deviation 
rate functions again obey the Hopf-Lax semigroup formula, as did the limit (|44p . 
The assumption is that lower tail rate functions exist initially: for all y, s G M 
the limit 

My,s) = - lim n-HogP{zJl m < ns} 



exists and is left continuous in y for each fixed s. Define 

^{w,r)^- lim ?i"MogP{r((0, 0), (n,nw)) < nr}. 



(52) 



This limit exists by superadditivity. For technical reasons a uniform tail bound 
is needed for the initial particle locations: there exist constants < Cj < oo 
such that 

P{zr(0) < -Ci|z| } < e-^^l'l for i < -Can, for large enough n. 
Theorem 4.2 ()Seppalainenl . Il998dl ) The limit 



Jt{x,r) = - lim n ^ logP{zJl^^Ant) < nr} 
exists for all x, r e R and t > 0, and is given by 

Mx,r) = ^ ^inf_Jjo(y,.)+t^(^,^)}. 



{y,s):y<x 



(53) 



24 



The approach of Section [5] can be adapted to prove the upper tail LDP 
((50|) and Theorem 14.21 The stationary systems make exphcit calculation again 
possible. Presently it is not clear how to include the lower tail LDP ([5T|) in 
the variational framework. Hence this part requires a separate proof. Details 



appear in Deuschel and Zeitouni ( 1999f ) and Seppalainen (1998d) 



EYirther remarks. Even though Theorem 14.11 has explicit rate functions, 
this large deviation problem remains unfinished in an important sense. It is not 
understood how the system behaves to create a deviation, and it is not clear 
what the rate functions / and U represent. The present proofs are too indirect. 

Let us illustrate through the random walk LDP (HI)-® how a large deviation 
problem ideally should be understood. To create a deviation Sn « nu with 
u > V, the entire walk behaves as a random walk with mean step u. Namely, 
it can be proved that the conditioned measure P( • | 5,i > nu) converges on the 
path space to the distribution Q'"^ of a mean u random walk. The value I{u) 
of the rate function in ^ is the entropy of this measure Q'^"^ relative to the 
original P. 



Re sults of the type presented in thi s section ap p ear fo r TA SEP in (ISeppalainenl . 
1998a ). The asymptotic analysis of iBaik et al. ( 19991 ) and Johanssonn 2000 ) 
gives also LDP's for the grow th models . Concentration results for a Brownian 



last-passage model appear in (jHamblv et al. and a general discus sion of 

deviation inequalities for growth models in the lectures of iLedouxl (|2007l) . 
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